Effects of nonlinear sweep in the Landau-Zener-Stueckelberg effect 
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We study the Landau-Zener-Stueckelberg (LZS) effect for a two-level system with a time- 
dependent nonlinear bias field (the sweep function) W(t). Our main concern is to investigate the 
influence of the nonlinearity of Wit) on the probability P to remain in the initial state. The dimen- 
sionless quantity e = 7rA 2 / (2hv) depends on the coupling A of both levels and on the sweep rate v. 
For fast sweep rates, i.e., e <C 1, and monotonic, analytic sweep functions linearizable in the vicinity 
of the resonance we find the transition probability 1 — P = e(l + a), where a > is the correction 
to the LSZ result due to the nonlinearity of the sweep. Further increase of the sweep rate with 
nonlinearity fixed brings the system into the nonlinear-sweep regime characterized by 1 — P = e 1 
with 7^1, depending on the type of sweep function. In case of slow sweep rates, i.e., e S> 1, 
an interesting interference phenomenon occurs. For analytic Wit) the probability P — Poe~ v is 
determined by the singularities of \/ A 2 + W 2 (t) in the upper complex plane of t. If W(t) is close to 
linear, there is only one singularity, that leads to the LZS result P — e~ E with important corrections 
to the exponent due to nonlinearity. However, for, e.g., W(t) oc t 3 there is a pair of singularities 
in the upper complex plane. Interference of their contributions leads to oscillations of the prefactor 
Po that depends on the sweep rate through e and turns to zero at some e. Measurements of the 
oscillation period and of the exponential factor would allow to determine A, independently. 

PACS numbers: 03.65.-w, 75.10.Jm 



I. INTRODUCTION 



We will study the Landau-Zener-Stueckelberg (LZS) 
problem of quantum-mechanical transitions between the 
levels of a two-level system at the avoided level cross- 
ing (see Fig n lL .caused by the time dependence of the 
HamiltonianEraHQEl 

H=~W(t)a z + ~Aa x , (1) 

where o~ a , a — x, y, z are the Pauli matrices and 

W(t)=E 1 (t)-E 2 (t) (2) 

is the time-dependent bias of the two bare (A = 0) en- 
ergy levels. Because of its generality, the LZS problem 
is pertinent to many areas of physics. In particular, the 
time-dependent model of Eq. (|]) is a simplification of the 
two-channel scattering problem with the curve crossing in 
the theory of molecular collisions (see, e.g., Refs. 0,0 and 
references therein). An appropriata-|iph©jce.,©f W(i), in- 
cluding oscillative time dependences,aiJ'tlJ'tilll3 may allow 
to m anipu late the evolution of the system in a controlled 
wayBOGj The most important case is probably that of 
the linear energy sweep: 

W{t) = vt, v = const > (3) 

that can be solved exactly. If the system at t — > — oo 
was in the bare ground state i/ji = | J.) before crossing 
the resonance, the probability to stay in this^istate after 
crossing the resonance at t — > oo is given by 



P(oo) = P = e 



ttA 2 
2hv 



(4) 
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FIG. 1: A pair of tunnel-splitted levels vs. energy bias 
W(t). Here E\ and E2 are the bare energy levels (A — 0), 
whereas E± are the exact adiabatic energy levels of Eq. (p3[). 
P denotes the probability to remain in the (bare) state 1 after 
crossing the resonance. 



matrix) and for multilevel systems can be found in Refs. 

EE 

Arguing that a typical smooth sweep function W[t) 
crossing the resonance-paUs one time can be linearized 
close to the resonancedOO and that the t-dependence 
of W{t) far enough from the resonance does not mat- 
ter, the LZS problem has been in most cases solved for 
linear sweeps Wit) = vt (see, however, Ref. |l8|). Re- 
cently, however, there has been a great deal of interest to 
the LZS effect in systems of interacting tunneling species 
(see, e.g._Etef. initiated by experiments on molecular 
magnets.E3 In Ref. [l9] the LZS effect has been investi- 
gated for an interacting system of N particles with spin 
S = 1/2. Since in this model each spin interacts with all 
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the others with the same coupling constant, one can use 
a mean- field approach which becomes exact for N — > oo. 
Therefore in general the effective field becomes nonlinear 
in time, due to the nonlinear i-dependence of the field on 
a given spin acting from other spins. 

It is thus of principal interest to investigate the two- 
level Landau-Zener-Stueckelberg problem for a nonlin- 
ear energy sweep. First, for a weakly nonlinear sweep 
there will be corrections to Eq. (||) that should be espe- 
cially pronounced in the slow-sweep region £>1. Since 
an experimental realization of a linear sweep can only 
be approximate it is important to determine these cor- 
rections. Second, essentially nonlinear sweeps, includ- 
ing those with multiple crossing the resonance, are be- 
coming technically feasible for such objects as magnetic 
molecules, where the Hamiltonian can be easily tuned by 
the external magnetic field (see, e.g., Ref. [H]). These 
nonlinear energy sweeps could be used to manipulate the 
system in a desired way, say to ensure P = for the finite 
sweep rates. 

In our recent LetterEl we have solved the inverse 
Landau-Zener-Stueckelberg problem, i.e., we have found 
such sweep functions W(t) that ensure the required time- 
dependent probabilities P(t) to stay in the initial state. 
Among the solutions are those corresponding to the time- 
symmetric P(t) satisfying P{— oo) = 1 and P(oo) = P = 
0. The functions Wit) are more complicated than cor- 
responding P(t) and they become cumbersome for more 
sophisticated forms of P(t). An interesting task is to find 
general criteria for the complete conversion tpi —> ip 2 , i-e., 
for P = 0, for the direct LZS problem and to apply it to 
the simplest nonlinear forms of W(t). This is a further 
aim of our paper. We will obtain analytical solutions for 
fast and slow sweep as well as numerical solutions in the 
whole range of parameters. 

The outline of the paper is as follows. In the next 
section we will describe the LZS problem for a general 
sweep. Sec. [II contains the discussion of the fast sweep, 
whereas results for slow sweep are presented in Sec. IV. 
We will show that in the slow-sweep case different kinds 
of singularities that characterize W(t) sensitively influ- 
ence the propability P. The final section contains a sum- 
mary and some conclusions. 

For numerical calculation throughout the paper we 
use Wolfram Mathematica 4.0 that employs a very accu- 
rate differential equation solver needed for dealing with 
strongly oscillating solutions over large time intervals. 



II. LZS PROBLEM FOR GENERAL SWEEP 

The Schrodinger equation for the coefficients of the 
wave function ip(t) = ci(t)tpi + c 2 (t)ip 2 can be written as 



lh Jt Cl{t) = T C2W ' 

ih J t ^(t) = -W(t)c 2 {t) + ^(t), (5) 



where c\ j2 {t) = Ci 2(i)exp 
the initial conditions 

ci (—00) = 1, 



2h 



C2 



0. 



satisfy 



(6) 



and W(t) = Ei(t) - E 2 (t) satisfies W(-oo) = -00. 
A general nonlinear sweep function W(t) can be 
parametrized as follows 

W(t) = Aw(u), u = vt/A. (7) 

The corresponding dimensionless form of Eq. (^|) is 

ie , 



c x {u) = -—c 2 (u), 

is 

c' 2 {u) = iew{u)c 2 (u) - — £i(u) 



(8) 



2 

— e. 

7T 



(9) 



where ' denotes differentiation with respect to u and the 
sweep-rate parameter e is defined by 

-= ^! 

hv 

(The reader should not confuse e and e that differ by a 
numerical factor; Whereas e arises in a natural way, e is 
used to represent most of the final results.) 

It is convenient to represent dimensionless analytical 
functions w(u) that behave linearly near u = in the 
form 

w{u) = a^fiau) *u+ |/> 2 + + . . . (10) 

where fk = 1 and, in general, /q ~ /q" ~ 1. Whereas v 
in Eq. ([]) stands for the sweep rate, the parameter a in 
( |l0| ) controls the nonlinearity of the sweep. Below we will 
work out explicit results for the analytical sweep function 



w(u) = a 1 sinh(au) 



11) 



that satisfies w(—u) — —w(u). A more general analytical 
sweep function is 

w(u) = — l — (e au - e~P u ) (12) 



that is characterized by w(—u) 7^ —w(u) for a 7^ j3. An 
interesting property of the corresponding LZS problem 
is that for a / (3 interchanging a and (3 leads to essen- 
tially different solutions of the corresponding Schrodinger 
equation while P being the same in both cases. This can 
be proven by considering the general scattering matrix 
of the problem, similarly to the proof of the reflection 
and transmission coefficients for scattering on a nonsym- 
metric one-dimensional potential baing the same for both 
directions of the ongoing particle.Ej 

Also we would like to consider nonanalytical sweep 
functions like the power-law w(u) = u a with a > 0. It is 
convenient to put this problem into a more general form 
and to introduce crossing and returning sweeps 



w(u) = 



sign(u)|w| a , crossing 
— I it I", returning. 



(13) 



Note that returning sweeps (with double crossing the res- 
onance) are naturally realized in atomic and molecular 
collisions where the distance between the colliding species 
at first decreases and than increases again. u 
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III. FAST SWEEP 

In this section we will discuss the behavior of P for 
fast sweep for which the propability P should stay close 
to 1. The form of Eq. (||) is convenient to perform the 
fast-sweep approximation e <C 1. In zeroth order of the 
perturbation theory one has ci(u) = 1 which can be used 
to obtain 62(11) from the second line of Eq. (||). The 
resulting equation for 62 (it) can easily be solved. Using 
this result, the probability to stay in the initial state 1 
can be expressed as 



pssl-l 

4 



5 2 (oo)| 2 

duexp 



-is I du'w(u) 
la 



(14) 



One can see that this is not a standard perturbation 
theory that would yield a correction of order e 2 to P. 
The integral over u in Eq. ( |l4| ) assumes large values 
for e <C 1 and is in general non-analytic in e. Since 
u ~ u e ~ £~ x / 2 1 contribute to this integral for e <C 1, 
the latter is very sensitive to the nonlinearity of w(u). 
For weakly-nonlinear w(u) that can be expanded into 
the Taylor series of Eq. (10) one has 



duw(u) = ^y - " 2 + 



3! 



4! 



If a <C u^ 1 ~ s 1 ' 2 , it is the first term of this expansion 
that is dominating. One can transform this condition 
into the dimensional form if one writes W(t) = Wot + 
Wot 2 /2\ + . . . and makes comparison with Eq. (|l0|). This 
yields the weak-nonlinearity condition 



a « s 1 ' 2 



nw Q 2 < wg. 



(15) 



With the help of Eqs. (g) and (|) one can establish that 
the LZS transition takes place not in the range W(t) = 
vt ~ A around the resonance, as could be expected, but 
in the much wider range 



W(t) ~ A/VI = Vfo > A 



(16) 



for the fast sweep. In the weakly nonlinear regime the 
main contribution to the integral in Eq. ( |l4| ) stems from 
u ~ u E and one obtains 



P^ 1 



1 + 



7T 

32 



f' 
Jo 



- f" 

, J 



111 



(17) 



Since a 2 <C e, the leading-order result for P does not 
depend on the nonlinearity, as expected. But Eq. Jl7j ) 
also demonstrates that the next-to-leading order in the 
nonlinearity reduces the probability P, independently on 
whether the sweep function w(u) [or W^(t)] grows slower 
or faster than linear, in contrast to the expectation. This 
effect that dependes on f' ' and /g" can be rather small, 
though (see Fig. ||). 

If a ^> e 1 / 2 and w(u) is a power series that terminates 
at finite order it™, then it is this latter term that dom- 
inates, and the main contribution to the integral in Eq. 
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FIG. 2: Dependence of the transition probability 1 — P on 
the nonlinearity parameter a for the sweep function w(u) = 
a -1 sinh(aii) in the fast-sweep regime, e < 1. The dashed 
lines are the asymptotes of Eqs. (JitI) and (E2J) . 



( |l4| ) comes from it e ~ e 1 /(™+ 1 ) that yields 



1 - P ~ (e/a 2 ) 



(18) 



The range of the energy bias W(t) that is responsible for 
tunneling is given by 



W(t) - A/s n/{n+1) > A 



(19) 



[cf. Eq. ©]. 

A particular case of the above is the crossing power- 
law sweep described by Eq. ( [To] ) that contains the linear 
sweep as a particular case. Here one has J Q duw(u) — 
M Q+1 /(a + 1) and Eq. (§3) yields 



p = 1 - r 2 



l + at Vl + a 



2a/(X+Qt) 



(20) 



where Y{x) is the gamma function. In the case a = 1 
the familiar expansion P = 1 — e of Eq. (|^) for e <C 1 
is recovered. For a > 1 the transition probability 1 — P 
is smaller than that for the linear sweep. The opposite 
result is obtained for a < 1, e.g., 1 — P oc \fe for a = 1/3. 



For the returning power-law sweep of Eq. (13) one has 
J duw(u) — — sign(w)|u| Q+1 /(a + 1) and the result has 
the form 



p = 1-r 2 



1 



2q/(1+q) 



l + aj~ \2(1 + a) J \1 + a / 

_ (21) 

For the sweep described by Eq. (O) one has 
/" duw(u) = (coshau — l)/a 2 - In the case a 2 <C e Eq. 
(O) yields the well known result for the linear sweep 
P = 1 — e, as it should be for small enough nonlinear- 
ity. In the opposite case one can use the approximation 
/ U du'w{u') e au /{2a 2 ) to obtain 



P S 1 - — ( In 



2a 2 



(22) 
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where 7 = 0.577216. This is in accord with Eq. (|18| ) in 
the limit n — > 00, up to the logarithm. The crossover 
between the linear and nonlinear regimes for the sweep 
function of Eq. (JTTI) on a/^/e is shown in Fig. 0. 



IV. SLOW SWEEP 

In this section we will study the asymptotic behavior 
of P for e » 1, i.e., for slow sweeps. This behavior 
will strongly depend on the analytical properties of w{u). 
Let us first describe the general situation in case of slow 
sweep. 



overbarrier-reflection problem with a complex potential. 
The result is a linear combination of two solutions that 
can be interpreted as ongoing and reflected waves. Then 
the (exponentially small) factor in front of the reflected- 
wave solution can be found from the analysis of Stokes 
lines in the complex plane (see, e.g., Ref. |^ and refer- 
ences therein). This analysis is rather involved, however. 
Below we will present an alternative and more simple 
method of solving Eq. ([?]) that does not rely on the WKB 
approximation. 

One can immediately write down the formal solution of 
Eq. (0) for |c + (oo)| 2 that contains yet to be determined 
c-(u) 



A. General 

It is convenient to solve the Schrodinger equation in 
the adiabatic basis formed of the states ip±(t) that are 
solutions of the stationary Schrodinger equation at the 
time t. The evolution of the system nearly follows the 
lowest, i.e., E- (t), of the adiabatic energy levels (see Fig. 

B> ' 

25±(t) = ±in(t) = ±\y/W*(t) + A*, (23) 

whereas the probability of the transition to the upper 
level E-u—is small. The explicit form of the adiabatic 
states isE3 

V±(t) - [±K ± (t)1> x + K T (t)^ 2 ] , 



where K±(t) = y/1 ± W(t)/fl(t). Writing the wave func- 
tion in the form if> = c+(t)ip+{t) + C-(t)ip-(i) one obtains 
the equations for c±(t) = cxp (i J dtE_{t)/ti) c±(t). The 
dimensionless form of these equations reads [see Eq. (]?])] 

c'_(u) — i - c+(u), fi(tt) = + w 2 (u) 
2fl (u) 

c'(u) = -ieU(u)£ + (u) - ^ 7 (M) c-(u) (24) 

2n (u) 

and the initial condition is c_(— 00) = 1, c+(— 00) = 0. 
The probability P to stay in state 1 is given by 



IMoo)| 2 
with 



- r (fu^^c_(u)exp [»£*(«)] (26) 
2 J -00 n (u) 



= / du'VLlu). (27) 
Jo 



Since for the slow sweep c_(u) remains close to 1, 
c_(ti) => 1 is a reasonable approximation. We will see 
below that this approximation is sufficient to obtain the 
correct exponential in the exponentially small P, as well 
as the prefactor with better than 10% accuracy. On the 
top of it, the approximation will be refined to obtain the 
exact prefactor. The asymptotic e dependence of the rhs 
of Eq. (^6|) depends strongly on the analytical proper- 
ties of $(u). Even if w(u) is chosen to be an analytical 
function, the integrand Q(u) in Eq. ( p7j ) is not analyt- 
ical but has singularities at the branch points at which 
w 2 (u) = — 1. We will see that these singularities deter- 
mine the large-e dependence of P, for analytical sweep 
functions w(u). 



B. Non-analytic sweep functions 

We consider for the beginning the simpler case of sweep 
functions that are non-analytic at crossing the resonance, 
u = 0, namely the power-law sweep functions of Eq. ( |T^ ) 
with a general a. For e ^> 1 the integral in Eq. (^6|) with 
C-(u) 1 is dominated by small u for which w(u) <C 
1, and thus O(tt) = 1. After taking advantage of the 
symmetry of Eq. ( |l3| ) it simplifies to 



p= f |2+(°o)| 2 , 

|c_(oo)| 2 = 1 — |c + (oo)| 2 , returning 



crossing , . 

2 \ lt> ) 



and it is small for crossing sweeps and close to 1 for re- 
turning sweeps. In fact, for returning sweeps P tends to 
1 in the limits of both fast and slow sweeps, as is illus- 
trated in Fig. ||. Indeed, from Fig. [l] one can see that for 
a fast sweep the system practically remains on the bare 
level E\ , whereas for a slow sweep it travels along the adi- 
abatic level P_ and thus returns to E\ for W(oo) = —00. 

One possible way of solving Eq. (Q) is to transform it 
to a single second-order differential equation and then 
to apply the WKB approximation to the corresponding 



|c+(oo)| =lim 
5^0 



due- Su au° 



i j cos (eit) 
sin(£w) 



for crossing and returning sweeps, respectively. The final 
result for P of Eq. J25|) reads 



P^ 



T 2 (l + a) o/"n-a\ 

— ¥°> — cos V~2/ ' crossm s 

1 r 2 (l + a) . 2 /?ra\ 

1 sin — , returning. 

e 2a V 2 / 



(28) 



The nonanalytic e dependent contribution here stems 
from the nonanalytic behavior of w(u) at u — and 
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FIG. 3: Staying probability P(e) for power-law sweep func- 
tions with exponents a — 1/2 and a = 1. Solid lines are 
numerical results and dashed lines are some of the analytical 
asymptotes. 



FIG. 5: -P(e) for crossing power-law sweep functions with 
exponents a = 1.5, 2 and 2.5. Dashed lines are high-e asymp- 
totes of Eq. Q. 




FIG. 4: -P(e) for power-law sweep functions with exponents 



it vanishes for analytic functions such as w(u) = u or 



w{u) 



for the crossing sweep and w{u) = u 2 for the 



returning sweep. Numerical results for P Tel along with 
the asymptotes of Eqs. (|2p|), (21), and (|28| ) are shown 
for a — 1/2 and a = 1 in Fig. H One can see that 
the slow-sweep asymptote of Eq. (|28| ) works well starting 
from relatively low values of e for crossing sweeps and 
from larger e for returning sweeps. P(e) for a = 1, 2, 3, 4 
and 5 are shown in Fig. In this case the asymptotes 
of Eq. (28) require rather large values of e and they can 
be seen on the log plot only (see Fig. ||) With increasing 
a oscillations of P(e) develop; the reason for this will be 
explained below. 



C. Analytic sweep functions 

As we will see, the form of the probability P for 
e » 1 crucially depends on the number of singularities 



of Q(u) = y/1 + w 2 (u) in the upper complex plane. In 
the next subsection we will discuss sweep functions w(u) 
for which there is only one singularity, starting from the 
weakly nonlinear case. For strongly nonlinear sweep with 
retardation in the resonance region, there is a pair of sin- 
gularities symmetric with respect to the imaginary axis, 
that dominate P for e 3> 1. In this case P can be repre- 
sented in the form of a small exponential with a prefactor 
oscillating as a function of e. 



1. One singularity: Exponent 

For analytic sweep functions w(u) the value of 
|c+(oo)| 2 given by Eq. ( |2^ ) is exponentially small for 
e > 1, and it is defined by the singularities of the in- 
tegrand in the upper complex half-plane that are closest 
to the real axis. These singularities are combinations of 

-1. Let us 



poles and branching points at w (u) 
consider at first sweep functions w(u) that are monotonic 
and close to the linear function, w(u) w u. In this case 
there is only one singularity at u ~ w{u) = w c = i, and 
it is convenient to use w as the sweep variable instead of 
u in Eq. ©: 



P^ 



2J-ooTf( W ) 



Jo 



dw 



dw 



(29) 



We deform the integration line into the contour C that 
goes down along the left side of the cut [i,zoo], then 
around the pole at w = i, and finally along the right side 
of the cut to ioo. Then the phase $(w) can be written in 
the form 



= $ c + 5$(w) 
du(w)- 



4> 



dw- 



dw 



fl(w) 
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6$(y) = i dyu'(y)y/V- y 2 , w 



ly. 



(30) 



Accordingly the result of integration in Eq. (g9|) can be 
represented in the form 



r- e 77 = 2elm$ c , 

where the prefactor Pq follows from 
1 f dy 



Po 



i - y' 



-c-{y)e 



(31) 



(32) 



Large values of e in our case ensure exponential smallncss 
of P in Eq. (|3l"|). Let us at first calculate the exponent 
r] for a weakly- nonlinear sweep described by Eq. (10). 
Inverting Eq. (fol) for a <C 1 one obtains 



One can notice that this prefactor slightlyjjtevid,tes from 
the exact prefactor 1 for the linear sweep.ofl This is the 
consequence of the approximation c_ (w) 1 in Eq. (|2S|). 

The reason for that wrong prefactor of Eq. ( |35| ) is the 
following. Although c_(w) only slightly deviates from 1 
along the real axis, it becomes singular at the relevant 
point w — i m the complex plane. This singularity can 
be worked out and the prefactor can be corrected. To this 
end, we note that the solution of Eq. ( pjj ) for e ^ 1 can be 
expanded in powers of 1 /e. This expansion captures only 
the analytical part of the solution that yields P = in all 
orders in 1 /e while the nonanalytic part of the solution 
that yields P ~ e~ e cannot be found by this method. 
Nevertheless, the 1/e expansion is sufficient to determine 
C-[w) that enters Eq. ( p2| ) to correct the prefactor. So 
we write down the expansions for c±(w) in the form 



of 



y 



c±{w) = 



C±, n (w) 



and accordingly 

rj = 2eTm$ c = si 1 



(36) 



(33) 



where c_ jn (— oo) = <5 n oi an( l c + .„(— oo) = 0. From Eq. 
(p3) that can be rewritten in the form 



dc-(w) c+(w) 



For the time-antisymmetric sweep, W(—t) — —W(t), one 
has /q = 0. Then for /g" > the exponent rj increases 
and the probability P to stay in state 1 decreases. This 
result seems counterintuitive since in this case the system 
spends less time in the vicinity of the resonance than in 
the case of a linear sweep and one could expect that P 
should increase. In fact, however, for a slow sweep the 
(small) probability to occupy state (+) oscillates many 
times during the resonance crossing and thus decreasing 
of P is a coherence effect that cannot be explaned by 
simple kinetic arguments. 



dw 

dc + (w) 
dw 



2n 2 ( w ) 



du(w) - 

-IE ; ll(w)C+(w) 

dw 



c_(w) 
2Tl 2 (w) 



(37) 



follows the infinite set of equations (u' = du{w)/dw) 



d5-, n (w) _ C +t n(w) 



dw 



2fi 2 
1 



dc +in -i(w) C- tV ,-x(w) 



dw 



2Q 2 



(38) 



2. One singularity: Prefactor 

Now we turn to the calculation of the prefactor Po m 
Eq. ( |3l| ) setting c_(y) => 1 in Eq. ( |32"| ) [sec comment 
below Eq. (p7|)]. For e^> 1 the values of y in the vicinity 
of 1 dominate the integrals. Thus one can introduce t = 
y - 1< 1 and simplify Eq. to 



This set of equations can be solved recurrently 
c +t o(w) = 0, c_ (w) = 1, 



1 



dw' 



,(39) 

2u'n" 4 J-oo u'Q," 

etc. In Figs. ||a and we compare the dependence 
|c+(w)| 2 obtained by the numerical solution of Eq. ( |37| ) 
and that using the n = 1 term of Eq. ( |36| ) with c+,i from 
Eq. ( |39"| ) for the linear sweep [w(u) = u] with e = 15. 
Whereas the agreement in the region w ~ can be im- 
proved by taking into account further terms of Eq. (pfl) , 
the corresponding smooth and even functions |c+(iy)|^ rt 

where u' c = u'(l) and the contour C t goes from oo to do * ot Jf^ ^ e asymptotic value P = |c+(oo)| 2 ^at^any 
above the cut [0, oo] around the pole at t = and then 
to oo below the cut. At the upper and lower sides of the 
cut one has (~t) 3 J 2 = ±it 3 / 2 . Thus Eq. @ yields 



Po 



dt 



■ exp 



(34) 



Po 



2?ri - 2i 



dt 

T 

2 



eu' c V2-t 



3/2 



1.096. 



(35) 



order of 1/e. The latter is due to the singular part of the 
solution that oscillates and tends to the exponentially 
small value P = e~ € [see Fig. 

To obtain the solution of the problem that is improved 
by use of the 1/e expansion while not losing the singular 
contribution at w — > oo, one should substitute Eq. ( [36] ) 
into Eq. (p9|). After multiple integration by parts the 
latter takes up the form 



Note that this prefactor is a number and it does not de- 
pend on the argument of sin, i.e., on e and on the con- 
stant u' c that encapsulates the nonlinearity of the sweep. 



P S 



n c-, n (w) 



71=0 



n 2 ( w ) 



(40) 
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FIG. 6: Lowest-order 1/e expansion for |c+(w)|p ert compared 
with the exact numerical solution. a-General view; 6-Right- 
wing magnification. 



I being the integration operator. The integration by 
parts eliminated the powers of 1/e that are present in 
Eq. (|36]). On the other hands, it is clear that contri- 
butions of the terms with n > in Eq. (40) should be 



small for e 3> 1 unless these terms have singularities at 
w = i. These singularitities do exist, as can be seen from 
Eq. ([39|) . So in the leading order in 1/s it is sufficient to 
take into account only the most singular contributions in 
c± n (w) near w = i that allows to simplify Eqs. ( |38| ) using 

w = i + z, \z\ <C 1, 

Q = V2iz, v! = u' c = const. (41) 

Then Eqs. ( |38| ) become 

dc- l7l (z) _ c +in (z) 



dz 



liz 

1 



uW2iz 



dc +t „-i(z) c_ i7l _i(z) 



dz 



4iz 



(.42) 



The solution of these equations can be searched in the 
form 



c±, n (z) 



P±,T 



(2i)"/2 („/)" z 3n/2 ' 



(43) 



where the coefficients j3±. n satisfy 

- 6n/3_,„ = P +t n, /?-o = 1 

0+,n = |(n-l)/3+, n _i + i/9_, n _i (44) 

and can be obtained recurrently. Under the same condi- 
tions, the sum in Eq. (J40|) simplifies to 



n— 

00 

-E 



n 

u'V2izI 



n=0 



2iz ^ V 3 

n=0 



(2i)«/ 2 +! «)" z ^/2+l 

i* 2^ I q J „i o,> 2^ d «' 



(45) 



n=0 



where we used Iz a — z a+1 /(a + 1) and defined 6 n = 
(2/3)™/?_,„/n! For S n one obtains from Eqs. Q) the re- 
currence relation 



5n = 1 - 



6n 



1 - 



S n -i, 



So = 1. 



Its solution is 



fe=l 



1 6fc 



(l/6)„(-l/6)„_(l/6)„(-l/6)„ 1 



where (a) n = T(n + a)/T(a). The sum X)^Lo ^« m ^1- 
(|45|) is a hypergeometric function of argument 1: 



E*« = E 



(l/6) n (-l/6) n »" 



(1), 



n 



n=0 n=0 

= a Fi(l/6, -1/6; l;x)| x=1 
1 3 

r(i + i/6)r(i- 1/6) ~ tt' 

Note that the approximation c_ (w) =>■ 1 in Eq. ( ^9| ) that 
leads to Eq. (|3^) for the prefactor amounts to neglection 
of all the terms of this sum except for the zeroth term 
5q = 1. Summing up all the most singular terms of the 
1/e expansion, as was done above, compensates for the 
wrong factor (7r/3) 2 in Eq. ( |35| ) and renders the prefactor 
the correct value 1. This value of the prefactor is valid 
for both linear and nonlinear sweep, if P is dominated by 
only one singularity in the complex plane. It is the con- 
sequence of cancelling of the quantity u' r , which describes 
the nonlinearity of the sweep, in Eqs. ( p5| ) and (fig). 

We thus have obtained the formula for the probability 
to stay in the initial unperturbed state 1 after a slow 
(e>1) energy sweep through the resonance 



P 



2s Im $ r 



dufl(u), 



(46) 
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where e is given by Eq. (||), D,(u) = y/\ + w 2 (u) is the 
dimensionless energy gap at the avoided level crossing. 
Note that 4> c does not depend on e. The integral in Eq. 
( |46| ) is taken from the real axis to the singularity point u c 
in the upper complex plane that is defined by Q(u c ) = 0. 
This formula can be found in the textbook by Landau and 
Lifshitz.113 Pokrovsky, Savvinykh, and UlinichcHI showed 
for a similar problem of the overbarrier reflection of a 
quantum-mechanical particle that the prefactor is equal 
to 1, using a different method. For a weakly- nonlinear 
swee p d escribed by Eq. ( JTo| ) the exponent n is given by 
Eq. (33). In contrast, as we have seen in Sec. [VB, for 



Im u 



non-analytic sweep functions W(t) the result is not expo- 
nentially small and is given by Eq. (^8|) for a particular 
case of crossing and returning power-law sweep functions 
ofEq. H). 



v\ w(u) = u 




1 1 1 — — 


,^j(3 m + 0/2 



Re u 



FIG. 7: Integration contour in the complex plane u in Eq. 

for the sweep function w{u) = it 3 . A pair of singularities 
at it = (±V3+*)/2 dominate the probability P for slow sweep, 
e > 1. 



3. A pair of singularities 

In many important cases the behavior of P is more 
complicated than the well-known formula Eq. (^) does 
suggest. For essentially nonlinear time- antisymmetric 
(crossing) or symmetric (returning) sweep functions W(t) 
typically there are two singularities at u c ± in the up- 
per complex plane that are closest to the real axis and 
symmetric with respect to the imaginary axis. It can 
be shown with the same method that P has the form 
(e»l) 



P 



$ c± = ±& c + i$'Z, (47) 



i.e., both contributions have prefactors 1. This is not 
surprising since for e 3> 1 contributions from different 
singularities are well separated from each other and the 
method used above applies to each of them. Eq. (47) can 
be rewritten in the form 



P 



P a e~\ P = 4cos 2 (eO 



n = 2i^ (48) 



1, 



that contains an oscillating prefactor Po with &' c 
according to the definition of the phase $ in Eq. (29) 
These results pertain to the crossing sweep. For the re- 
turning sweep one obtains 



P = 1 - P e- 



(49) 



Turning the prefactor Po to zero for some (finite!) 
sweep rates is the so-called complete conversion from 
state 1 to state 2 that was mentioned in the introduc- 
tion. In Ref. [l^ we have found this phenomenon for spe- 
cially chosen sweep functions W(t) of a more complicated 
form. Here we demonstrated that the full conversion at 
finite sweep rates is a general phenomenon for a nonlin- 
ear sweep. Whereas an oscillating prefactor can be found 
in earlier publications (see, e.g., Ref. |l^), here we have 
shown its simple relation to the relevant singularities of 
the LZS problem. 



D. Particular cases 



Let us now consider particular cases of the LZS ef- 
fect with nonlinear energy sweep, which will allow to 
explore the role of the singularities in the prefactor in 
more detail. For the power-law sweeps of Eq. ( |l3| ) with 
natural a — n, such as w(u) = u 2 , w(u) — vr, etc., 
there are n singularities in the upper complex plane at 
u = iik = exp [z7r(l/2 + k)/n] with k = 0, 1, .. .,n — 1 
(see Fig. for a — 3). Closest to the real axis are those 
with k = and k = n — 1. The integration contour in 
Fig. ^ goes along both sides of the cuts initiating at the 
singularity points. These cuts are chosen so that (5$ in 
Eq. ( j30| ) is real on both sides of the cut (i.e., the cuts are 
the so-called anti-Stokes lines), so that P is determined 
by purely oscillating integrals. For all other directions of 
the cuts, one has exponentially converging integrals on 
one side of the cut but exponentially diverging integrals 
on another side of the cut. After calculation of <& c in Eq. 
(Effl one obtains 



2n 

7T 

2n 



,1/n 



r 2 (i + M 



T 2 



(50) 



Note that for n = 1 there is only one pole, and Eq. (|50| ) 
yields <$>' c = and n = 2e$^' = (vr/2)e = e, so that 
one returns to the LZS result P = e~ e . In other cases 
the prefactor Po like in Eqs. (|4|) and ([l9]) oscillates and 
turns to zero at some values of the sweep rate, where 
contributions of the two singularities in Eq. fl47| ) cancel 
each other. In the latter complete Landau-Zener- 

Stueckelberg transition is achieved for crossing sweeps. 
We have shown the prefactor Po for w = u 3 in Fig. ||, 
where the solid line is e^P with numerically computed P 
and n = 2e$" and the dashed line is the analytical result 
Po = 4 cos 2 (e$;). It is seen that Eqs. © and © work 
well for large e. Plotting Pn reveals oscillations that are 
not seen at large e in Fig. [|. 

The next interesting case is that of the sinh sweep func- 
tion of Eq. ( pd| ) . This model allows to track the changes 
of P as function of the continuous parameter a (which 
is the measure of nonlinearity) keeping the sweep func- 
tion w(u) analytical. The singularities depend on a. For 
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FIG. 8: Prefactor Po in P{s) = Poe~ v for power-law sweep 
function w — u z . The solid line is e^P with numerically com- 
puted P and rj = 2e$" an d the dashed line is the analytical 
result P = 4 cos 2 (e$c) [see Eqs. (H and 



w(u) = a sinh(aw) 




FIG. 9: Real and imaginary parts of the quantum-mechanical 
phase $ at symmetric singularities in the complex plane [see 
Eq. (^)] for the sweep function w = a -1 sinh(au) vs a, given 
by Eqs. @ and (H). 



(51) 



a < 1 the singularity closest to the real axis is at 

arcsincv 



u = u r = i 



whereas the next singularity is at u = i (n — arcsina) /a. 
Keeping only the singularity at u c for e 3> 1 one obtains 
from Eq. (46) 



2e ( 1 
rj = — E ( arcsina, — ) , (a < 1), (52) 



where E(<y9, m) is the incomplete elliptic integral of second 
kind. In the region a > 1 there is a pair of singularities 
at 



u = u c ± = ± 



a \ 



a 2 — 1 + a 



in 



(53) 



a \ /2a 
that are closest to the real axis. In this case one obtains 



FIG. 10: Prefactor Po in P — Poe~ v for the sweep function 
w = a" 1 sinh(aii) defined as Po = e v P with numerically com- 
puted P for a in the vicinity of 1. The dotted lines at Po = 1 
and at Po = 2 correspond to the analytical results for a < 1 
and a = 1, respectively. 
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w(u) = (e -e p )/(a+p) 




FIG. 11: Staying probability P vs the inverse sweep rate 
for the sweep with slightly broken antisymmetry w(—u) = 
—w(u). For a = 9 and /3 = 11 (same results for a = 11 and 
[3 = 9) interference of the two poles closest to the real axis is 
not complete and P does not turn to zero. 



Eq. 



j) with 
V = 
& = 



2 ^ / 
— E 


a 2 


Q 




1 


K 




a 







E 1 



(54) 



Note that &' c = for < a < 1. The a dependences of 
77 and above are shown in Fig. ||. In the special case 
a = 1 both singularities join in one and the phase <!>(it) 
of Eq. (|2^) becomes a holomorphic function of a complex 
argument: 

$(u) = w(u). 

In this case one could be puzzledJai-the attempt to apply 
the general Landau's argumentsBEZl leading to Eq. (J4q) . 
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Nevertheless, the pole in Eq. ( |26| ) due to I/O. (it) remains 
and it determines the exponentially small P. It is again 
convenient to consider w as the sweep variable and to use 
Eq. (p9[). The coefficient c_(u>) near the pole at w = i 
is obtained from Eqs. (|37|) with the Ansatz of Eq. (36), 



this time taking into account the singularity of du/dw = 
1 / Vl + w 2 . As the result one arrives at Eq. ([l8|) with the 
prefactor 



Po = 2, (a = 1). 



(55) 



Oscillations of the prefactor Po m the region a > 1 
for this model can be illustrated in the same way as was 
done in Fig. [|. More interesting here is to plot P in 
the vicinity of a = 1 where the crossover between oscil- 
lating and non-oscillating regimes takes place (see Fig. 
0). Although the asymptotic numerical behavior of Po 
for a < 1 and a = 1 at e 1 is in accord with the 
analytical results above, one can see deviations from this 
picture at smaller e because of the contribution of the 
more distant pole at a < 1. 

One can also consider sweep functions that are not 
time-symmetric or antisymmetric, such as Eq. (12). If 
a and (3 slightly deviate from each other, there are two 
singularities that are at slightly different distances from 
the real axis. In this case one can observe oscillations 
of P in some range of e without exactly turning to zero 
(see Fig. [n]). At larger e the singularity closest to the 
real axis dominates and oscillations disappear. Note that 
the probability P is symmetric with respect to the inter- 
change of a and f3 in Eq. P2). 



V. CONCLUSIONS AND SUMMARY 

Our main concern has been the investigation of the 
LZS effect for nonlinear sweep. There were three main 
points we wanted to clarify. First, what are the correc- 
tions to the LZS propability P, due to weak nonlinearity 
of analytical functions W(t)7 This question is of primary 
interest since an experimental sweep can be linear only 
approximately. Second, are there any qualitatively new 
features for analytic but strongly nonlinear W(t) and, 
third, how does P look like for nonanalytical sweep func- 
tions? Of course, the answers to these questions depend 
on whether the sweep is fast or slow. 

Concerning the first point, we have found for fast sweep 
rates, i.e., e -C 1, that e sets the scale for the rele- 
vance of nonlinearities of W(t) measured by a dimen- 
sionless parameter a of Eq. Jl0|). For fast sweep and 
weak nonlinearity a 2 <C e <C 1, one obtains the ex- 
pected result 1 — P = e with positive corrections of order 
(a 2 /e) . For still faster sweeps and/or stronger nonlin- 
earities, e <C min(l,Q! 2 ), and for w(u) that is a power 
series terminating at finite order u n , the leading-order 

result is 1 — P ~ (e/a 2 ) 2 ™^™ +1 ' ) . This implies that the e- 
dependent transition probability 1 — P for the fast sweep 
changes from e for a linear sweep to s 2 for a nonlinear 
sweep with n — > oo, in leading order in e. This conclu- 
sion, of course, also follows from the result of Eq. (|2C 



for the pure power-law sweep defined by Eq. ( |l3|) since a 
corresponds to n. 

The most interesting results follow for slow sweeps, i.e., 
e 3> 1. For several qualitatively different sweep func- 
tions W(t) we have shown that P can exhibit quite differ- 
ent e-dependence. For the power-law sweep of Eq. ( |T3] ) 
which is a nonanalytical function one obtains a power- 
law behavior for P, i.e., P ~ e~ 2a in case of the crossing 
sweep. It would be interesting to experimentally real- 
ize such a sweep because it should be much easier to 
measure the e~ 2a dependence than an exponential de- 
pendence e~ e , for large e. For analytical sweep functions 
we have demonstrated that P can always be decomposed 
into an exponential factor e~ v = exp (— 2elm $ c ) and a 
prefactor Po = 4 cos 2 (eRe^c) , where $ c is determined 
by the singularities u c following from 



w 2 (u c ) 



(56) 



If w(u) does not deviate strongly from the linear func- 
tion [see the exact criteria in the main text] then there is 
one singularity in the upper complex plane, only, and one 
obtains the LZS result with Po = 1 and rj = e. However, 
a new feature occurs if the energy sweep W(t) exhibits 
a significant retardation in the vicinity of the resonance. 
In this case we have shown that the prefactor Po becomes 
oscillatory as function of s and turns to zero at some e. 
The latter corresponds to the complete conversion from 
state 1 to state 2ll3 Since the period of the oscillations 
of Po is proportional to e -1 ~ v/A 2 , its measurements 
would allow independent determination of A for a given 
sweep rate v, besides the measurement of the exponential 
factor. Of course, to determine A from the oscillation's 
period one must calculate Re $ c which is determined by 
W(t), only. We mention that oscillations of P were al- 
ready found for different nonlinear-sweep models (see, 
e.g., Ref. |l^), as well as for the exactly solvable tight- 
binding electron model on one-dimensional chain driven 
by an electric field .o 

Oscillations of the prefactor Po in Eq. ( p8| ) for es- 
sentially nonlinear energy sweeps find physical expla- 
nation. Consider, for instance, the sinh sweep of Eq. 
(|i"T|). The system is in the vicinity of the resonance for 
— uq < u < Mo, where uq satisfies w(uq) — 1 and is given 
by uq — In (Va 2 + 1 + a) /a. For a >• 1 one has u <C 1 

and the derivative w'(uq) — cosh(aito) = Vex 2 + 1 is 
large, thus w(u) < 1 in the main part of the interval 
— uq < u < uo- This means that for a 3> 1 the system 
rapidly comes into the vicinity of the resonance, stays 
practically at resonance for some time, and then rapidly 
leaves the resonance region. During this stay, the system 
can oscillate between states 1 and 2, so that the value of 
P depends on the time spent at resonance and at some 
values of e it turns to zero. The time spent at resonance 
to = (A/u) uq [see Eq. (Q)] is controlled by the sweep rate 
v or by the parameter e = A 2 /(h~v). Note, however, that 
for finite a oscillations between states 1 and 2 are not 
complete and thus P < 1 and it becomes exponentially 
small for e>1. The same physical explanation is valid 
for power-law sweeps w(u) — u n with n 1, although 
the effect is already seen for n > 1. 
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Similar situation should be realized in the case of the 
overbarrier reflection of a quantum-mechanical particle. 
If the potential barrier is parabolic (which corresponds to 
the linear sweep in the LZS problem) then the problem 
has an exact solution and the reflection coefficient R de- 
termined by a single singularity in the complex-plane is a 
monotonic function of the particle's energy Ex-H If, how- 
ever, the potential U(x) has a flat top, say U(x) = —x 2n 
with n = 2, 3, ... , then there is a pair of singularities 
nearest to the real axis that causes R(E) to oscillate and 
turn to zero for some values of E. This behavior is well- 
known for the rectangular potential that is the limiting 



case of U(x) = —x 2n with n — > oo. 

Finally, the 1 Je expansion of Eq. (|3^) that was applied 
to obtain the prefactor in the LZS formula for analyti- 
cal sweep functions and slow sweeps might be interesting 
from the technical point of view. 
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